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Objectives 

The  objective  of  this  research  is  to  develop  a  theory  for  the  behavior  of  a  rotor  wake  undergoing 
dynamic  perturbations.  The  theory  should  be  capable  of  finding  the  flow  in  all  three  directions  both  on 
the  rotor  disk  and  off  of  the  rotor  disk.  The  theory  should  be  in  a  state-space  form  that  is  hierarchical 
and  can  be  used  in  real-time  simulation  studies.  It  should  include  dynamic  inflow  and  the  He  Dynamic 
Wake  Model  as  special  cases,  and  should  be  derived  from  first  principles.  It  should  allow  for  mass 
source  terms  as  well  as  pressure  differentials. 


Approach 

The  approach  taken  is  as  follows: 

1. )  Formulate  equations  based  on  a  Galerkin  approach  on  the  linearized  potential  flow  equations 

for  a  rotor  in  a  skewed  free-stream. 

2. )  Develop  all  Galerkin  integrals  in  closed  form  based  on  the  divergence  theorem. 

3. )  Study  the  convergence  of  the  model  for  test  cases  with  known  solutions. 

4. )  Develop  strategies  to  overcome  singularities. 


Accomplishments 

The  history  of  dynamic  wake  modeling 
covers  the  last  55  years.  References  [1] 
through  [17]  outline  the  work  up  to  the  1990’s. 
Several  important  developments  followed  that 
changed  the  course  of  dynamic  wake  modeling. 
First,  there  was  an  unusual  anomaly  in 
helicopter  flight  mechanics  that,  during  a 
pitching  maneuver,  experimental  flight-test  data 
showed  that  rotorcraft  rolled  opposite  to  what 


the  simulations  predicted  even  with  the  new 
wake  models.  Rosen,  Ref.  [18],  suggested  that 
this  might  be  due  to  the  effect  of  wake 
curvature  on  induced  flow  during  the 
maneuver. 

Curtiss  and  coworkers  investigated  this 
and  found  that  this  curvature  effect  could 
explain  part  of  the  anomaly  and  that  simple 
vortex  and  momentum  considerations  could 
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capture  the  effect,  Refs.  [19-20].  Barocela, 
corroborated  these  findings  in  his  work  on 
wake  distortion,  Ref.  [21],  although  some 
discrepancies  were  noted  with  the  work  of 
Curtiss,  Ref.  [22].  Barocela  showed  that  the 
effect  of  wake  curvature  could  be  incorporated 
into  the  dynamic  wake  model  of  Pitt  and  Peters 
by  allowing  the  coefficient  matrices  to  be 
functions  of  both  the  wake  skew  and  the  wake 
curvature.  This  was  completed  by  Krothapalli, 
et.  al  in  Ref.  [23]. 

The  work  on  wake  curvature  was 
extended  to  the  He  model  in  References  [24-26] 
in  which  it  was  shown  that  the  entire  coefficient 
matrix  (involving  all  harmonics)  could  be 
modified  to  account  for  wake  curvature. 
Reference  [27]  demonstrated  that  the  work  of 
Curtiss  could  be  considered  a  special  case  of 
this  more  general  formulation,  thus  uniting  all 
of  the  various  formulations  together. 

Another  area  of  work  on  dynamic  wake 
came  with  the  need  for  ground  effect 
calculations.  The  original  He  theory  had 
included  a  simple  ground-effect  correction,  but 
it  became  clear  that  something  more  general 
was  needed  in  order  to  accommodate  moving 
ground  planes,  inclined  ground  planes,  and 
partial  ground  planes.  References  [28-31] 
progressed  through  a  series  of  developments  in 
which  the  He  inflow  model  was  used  to  treat 
ground  effect.  In  these  various  approaches, 
either  an  image  rotor  or  a  ground  plane  source 
rotor  was  used  to  simulate  the  ground  in  a 
quasi-steady  manner.  For  these  cases,  a  second 
actuator  disk  (for  the  image  rotor  or  ground 
rotor)  was  used,  and  the  induced  flow  from  that 
rotor  needed  to  be  computed  on  the  primary 
rotor.  This  was  done  numerically  by  an  off-line 
quasi-steady  approach. 

From  that  work,  it  became  clear  that  what 
was  needed  was  a  completely  unsteady  theory 
that  could  compute  the  induced  flow  of  these 
secondary  rotors  above  the  disk  plane  and  in  all 
three  components.  The  He  model,  however, 
was  unable  to  do  this;  and,  therefore,  the  need 
was  established  for  a  more  general  treatment  of 
inflow  dynamics.  This  led  to  an  entirely  new 
methodology  for  induced  flow  dynamics. 


Reference  [32]  outlined  an  approach  for 
deriving  dynamic  wake  equations  in  terms  of 
velocity  potentials  rather  than  merely  normal 
velocity  components  at  the  rotor  disk.  Such  a 
formulation  was  shown  to  be  derivable  from  a 
Galerkin  procedure  applied  to  the  potential  flow 
equations.  In  Reference  [33],  this  methodology 
was  successfully  implemented.  Reference  [34] 
demonstrated  that  this  method  contained  the  old 
He  model  (and,  consequently,  the  Pitt-Peters 
model)  as  a  special  case  in  which  off-disk 
modes  were  neglected.  References  [35-36] 
demonstrated  that  the  method  did,  indeed 
converge  for  all  three  components  of  flow  in 
the  entire  half  space  including  the  plane  of  the 
rotor  disk  and  the  semi-sphere  above  it. 

Fluid-Dynamics  Equations 

The  three-dimensional  potential  flow 
equations  (momentum  and  continuity 
equations)  for  the  pressure  and  velocity  fields  P 
and  v ,  with  a  free-stream  velocity  Vcc,  are 

f :<■> 

V»  v  =  <?  (2) 

These  equations  have  been  non- 
dimensionalized  by  defining  P  as  pressure 
divided  by  pVj,  v  as  induced  velocity  divided 
by  Voo,  and  time  as  a  reduced  time  r,  (i.e.,  time 
multiplied  by  VJR.)  The  variable  £is  the  non- 
dimensional  coordinate  along  the  free-stream 
line,  positive  upstream.  All  lengths  are  divided 
by  the  rotor  radius  R.  Figure  1  shows  the 
coordinate  system. 

From  continuity,  Eq.  (2),  it  is  observed 
that  v  can  be  expressed  by  a  velocity  potential 
that  satisfies  Laplace's  equation.  It  can  also  be 
shown  that  P  satisfies  Laplace’s  equation. 
Therefore,  P  can  be  expressed  as  a  summation 
of  pressure  potentials,  0 ;  and  v  can  be 
expressed  as  a  summation  of  the  gradient  of 
velocity  potentials,  lF. 
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Pressure  Potentials,  and  Velocity 
Potentials 

To  transform  Eqs.  (1)  and  (2)  by  a 
Galerkin  method,  it  is  required  to  expand  the 
pressure  potential,  <P ,  and  the  velocity 
potentials,  W ,  in  terms  of  a  complete  set  of 
functions,  each  of  which  satisfies  Laplace's 
equation.  In  addition,  they  have  to  fulfill  the 
boundary  conditions  for  pressure  in  the  case  of 
0 ,  and  for  velocity  in  the  case  of  ¥ . 

The  boundary  conditions  for  pressure  are 
given  by  a  discontinuity  across  the  rotor  disk, 
Fig.  (1).  The  use  of  an  ellipsoidal  coordinate 
system  [ v,7],y/  ],  has  the  advantage  that  any 
odd  function  in  v ,  will  allow  a  representation 
of  a  discontinuity  across  the  rotor  disk.  An 
additional  advantage  of  using  an  ellipsoidal 
coordinate  system  is  that  an  analytical  solution 
of  Laplace’s  equation  is  known  and  can  be 
expressed  as, 

(v,  T],  w)  =  P„m  (v}Qnm  (j  v)C0S{mW)  (3) 

Or  (v, T],  y/)  =  P„m {v)Qn  (i^)sin(m^)  (4) 

where  P„m(v)  and  Qn(irj)  are  normalized 
associated  Legendre  functions  of  first  and 
second  kind. 

Since  P„(v)  with  n+m  odd  is  an  odd 
function  of  v  as  well  as  a  function  that  satisfies 
Laplace's  equation,  it  has  a  discontinuity  across 
the  disk.  Therefore,  Eqs.  (3)  and  (4)  can  be 
used  as  the  expansion  functions  for  the  pressure 
potentials.  On  the  other  hand,  P„m  (v)  with  n+m 
even  is  an  even  function  of  v ,  but  its  derivative 
with  respect  to  z  is  an  odd  function  of  v  on  the 
disk.  Such  functions  can  be  used  to  represent 
mass  sources  at  the  rotor  disk. 

P  =  -£  z(c<c+C<“)  (5) 

m=0n=m+l 

The  boundary  conditions  for  the  velocity 
field  are:  (1)  the  velocity  field  far  upstream 
from  the  rotor  is  equal  to  zero,  and  (2)  there  is  a 
velocity  discontinuity  any  place  a  vortex  or 
vortex  sheet  exists  in  the  flow  field.  These 
discontinuities  only  exist  at  the  rotor  blades  and 
within  the  rotor  wake. 


If  the  velocity-field  computational 
domain  is  limited  to  the  infinite  upper-half 
volume  above  the  rotor  disk  (Fig.2)  the 
functions  to  be  used  in  the  velocity  expansion 
do  not  have  to  fulfill  any  discontinuity 
conditions.  It  is  important  to  note  that,  in  the 
case  of  perfectly  edgewise  flow,  the  wake  is 
located  on  the  rotor  disk  plane.  Thus  no 
convergence  of  this  methodology  is  expected 
on  the  trailing  region  off  the  rotor  disk  for 
edgewise  flow,  since  the  assumption  for  the 
velocity  potentials  is  no  longer  valid  in  this 
region. 

In  order  to  strongly  satisfy  the  upstream 
boundary  condition  for  the  velocity  field,  the 
velocity  potentials  are  defined  as 

kc = ;  r? = 7<w 

f  i 

m  =  0, 1,2,. n  =  m,m  +  l,m  +  2,...,°°  (6) 


This  strongly  ensures  a  zero  velocity  upstream. 
As  £  approaches  infinity,  ¥nm  approaches  zero. 
The  velocity  v  is  taken  as 

v=  £  £ +  ffi'F™)  (7) 

m-On-m 


where  ¥™c  and  are  defined  by  Eq.  (6). 

Substitution  of  Eqs.  (5)  and  (7)  into  (1), 
yields  the  momentum  equation  in  an  expansion 
form.  Because  a„  ,  b ”,  t™c and  r^are 
functions  of  nondimensional  time  (f)\  while 
ipnmc ,  ,  &™c ,  and  are  functions  of  the 

spatial  coordinates  ( x ,  y,  z),  Eq.  (1)  takes  the 
following  form  for  cos  ( )  terms: 


oo  oo 

I  I 

m-0 n-m 


-  da™  ~d 

y  tj/ttic  _  y  n 


dr 


*m 

-a„ 


(8) 


I  £v<crr 

m-0  n~m+l 

Only  the  cosine  terms  are  listed,  since  sine  and 
cosine  completely  decouple.  However,  an 
identical  set  of  equations  can  be  written  for  the 
sine  terms. 


Galerkin  Approach 

If  a  Galerkin  approach  is  applied  to  solve 
the  conservation  of  momentum  equation,  Eq. 
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(8)  is  pre-multiplied  by  the  gradient  of  each  one 
of  some  test  functions,  A,  integrated  over  the 
domain;  and  the  integrals  are  set  to  zero.  In  a 
Galerkin  methodology,  the  test  functions  are 
defined  from  the  same  set  of  functions  used  to 
expand  the  pressure  potential,  <P ,  and  the 
velocity  potentials,  5P.  Thus,  we  choose 

Aj  =<Prf\  Aj  =<Pjs 

r  =  0,1,2,...,°°;  j  =  r,r+  l,r  +  2,...,°°  (9) 

When  this  Galerkin  approach  is  applied,  the 
conservation  of  momentum  equation,  Eq.  (8), 
becomes 


JJT 


( 


■ZZ 

m-0n=m 


■d%_ 

dr 


V7  d\n_gn 


ty  = 


V  m=0n=m+l 

r  =  0,1,2, j  =  r,r  +  l,r  + 


(10) 


The  Divergence  Theorem  allows  volume 
integrals  to  be  expressed  as  surface  integrals, 
Fig.  2.  Therefore,  if  this  theorem  is  applied  to 
Eq.  (10),  and  because  Aj,  ,  and  all 
fulfill  Laplace’s  equation,  Eq.  (10)  becomes 


oo  oo 


,rcWZ  da.:  8* 

Aj  1 r~dT  j  dndi  n 


EZJI 

m=0 n-rn  $  V 

°o  oo 

Z  Z 

m=0n=m+l  S 


IS  = 


r  =  0,1,2, . . . ,  oo;  j  =  r ,  r  + 1,  r  +  2, . . . ,  oo 

or 

OO  oo  (  3\rf 

V  V  ff  ™mc  C'Aj  <^n  an 
00  00  dAr? 

I  £  JM-  KcdsC 

m=0n=m+l  5  </n 


(ID 


r  =  0,1,2,...,°°;  j  =  r,r  +  l,r  +  2,...,° 


(12) 


where  S  represents  the  surface  area  of  the 
upper-half  inflow  volume  V,  including  the  rotor 
disk  plane,  dS  represents  a  differential  area 
element;  and  n  is  a  unit- vector,  outward  normal 
to  dS.  The  surface  S  can  be  subdivided  into 
three  areas:  two  of  them  are  located  on  the  rotor 
disk  plane  (z=0),  s,  and  correspond  to  the  on- 
disk  area,  A,  and  the  off-disk  area,  B\  and  the 
third  one  corresponds  to  the  area  on  the  infinite 
dome,  C  (Fig.2).  Therefore, 


s  =  A+B;  S  =  s+C  (13) 

The  pressure  potentials,  <P™C ,  velocity 


potentials,  and  test  functions,  Aj ,  are 

such  that  all  the  integrals  over  the  dome 
surface,  C,  become  zero.  Therefore,  the  surface 
of  integration  becomes  the  rotor  disk  plane 
( z=0 ),  s,  and  the  normal  outward  vector  is 
along  the  z  axis.  Therefore, 


d_=d_ 
dn  dz 


(14) 


If  Eq.  (14)  is  substituted  into  Eqs.  (11)  and 
(12),  together  with  the  definition  of  the  velocity 
potentials  and  the  test  functions,  Eqs.  (6)  and 
(9),  they  become 


m=0n=m  s 


Z  £  K  Tz 


"n 


dz  dz 


oo  oo 

Z  I 

m=0n-m+l  $ 


r  =  0,1,2,...,°°;  j  =  r,r+l,r  +  2,...,°° 

f  °°  j *m  | 

“  „  uun  ^  mr  Am  » 

ds  = 


oo  oo 


,34'' 


Z  Z  JHf  S 

m=0n~m  s  °Z  \0 

oo  oo 

Z  Z 

m-0n=m+\  s  4 


(15) 


r  =  0,1,2,...,°°;  j  =  r,r  +  l,r  +  2,...,°° 


(16) 


where  ds  is  a  differential  area  on  z=0. 

At  this  point,  it  is  important  to  note  that  the 

with  n+m  odd  are  zero  on  the  z-0  plane  off 
the  disk  (B),  whereas  the  z  derivatives  of 
with  n+m  even  are  zero  on  region  B.  Because 
of  this,  the  Divergence  theorem  can  be  utilized 
with  appropriate  choice  of  the  3 /dz  position  in 
each  individual  term  such  that  all  integrals  are 
zero  on  region  B.  The  result  is  a  set  of  integrals 
that  need  only  be  evaluated  on  region  A  (on- 
disk)  for  which  they  can  be  evaluated  in  closed 
form.  The  expression  obtained  after  applying 
this  procedure  can  be  condensed  in  the 
following  equation. 

[ff]{ar}+[D']{a:}=[D']{c}  <17> 

where  (  )*  =  d(  )/dt  ond  where  each  one  of  the 
elements  of  the  lc  and  Dc  matrices  are  known 
in  closed  form. 
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Equation  (17)  is  valid  for  any  skew  angle 
which  appears  in  the  equation  in  the 
expressions  for  the  wake  influence  coefficient 
matrix  Lc .  This  equation  can  be  further 
partitioned  into  two  row-groups  and  two 
column-groups  such  that  m+n  (or  j+r)  is  odd 
and  m+n  (or  j+r)  is  even.  These  matrices  are 
organized  in  the  following  way 


If  Eq.  (17)  is  organized  as  suggested  in  Eq. 
(18),  it  can  be  partitioned  as 


[<X,o  !  [^L,e  jp“L  \p\),0  j  \p\>{ 

h;tec;J{^  } + lk7tm;; 

rofiRTlkll 


(19) 


Potential  Function  Expansions 


The  non-dimensional  pressure  drop  and 
mass  flow  added  to  the  velocity  field  (both 
across  the  disk),  and  the  velocity  everywhere  in 
the  upper-half  plane  can  be  computed  for  the 
cosine  terms  as 


py2  ~  \_Plower  Pupper  \T]=0 

2Z  Z  fniV)  [1T)0COS(”^) 

m~0  n-m+l,m+3,.... 


(20) 


—  = 
pVc o 


\ower  ***  Supper 


\=0 


2  Z  i  Pnm{ArnC)eCOSM 

m=0  n-m+2,m+4,.... 


(21) 


v=i  l  a"viCc=  i  £  far#  (22) 

m-0  n=m  m=0n=m  ^  j 

From  Eq.  (22),  it  is  seen  that,  to  compute 
the  velocity  field,  it  is  required  to  compute  the 
velocity  potentials,  ,  by  a  numerical 
integration 


K  =  \®nCdZ  (23) 

£ 

To  avoid  numerical  integration  and  in  order 
to  be  able  to  express  the  velocity  potentials  in 
terms  of  potentials  known  everywhere  in  the 
flow  field,  a  change  of  variable  from  a™  to  a” 
is  introduced 

k  F  k“l=  W  F  + 1  (24) 

The  constants  o™  and  g™  are  chosen  such  that 
the  new  velocity  potentials  will  give  no 
singularities  when  gradients  are  taken. 

(25) 


and 

[n)H  =  (n)(n-2)(n-4)...(2),n  =  even 

(«)//  =  {n)(n-2)(n-4)...(l),n=odd 

(0)!!  =  1;  {-1)U  =  1;  (-2)//  =  °°;  {-3)!!  =  -l 


(27) 

(28) 

(29) 


If  a  Galerkin  approach  is  applied  to  Eq. 
(24),  the  following  relationship  can  be  obtained 

fc!4‘FVk)  <30> 

where 

H=|ZCU  (3D 

and  where  |lc]  is  the  same  matrix  as  defined  in 
Eq.  (17)  and  in  the  Appendix. 

Equation  (24)  can  be  used  to  express  the 
velocity  potentials  in  axial  flow.  [32]  If  Eq. 
(24)  is  substituted  into  Eq.  (22),  it  yields 

v;=  I  i^vk^+srr^s)  (32) 

m-0  n=m+l 


the  axial  component  of  the  velocity  is  given  by 
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VZ=2  I  (34) 

m-0  n-m+1 

We  recall  that  the  He  model  came  in 
two  equivalent  versions,  one  with  velocities 
expressed  in  terms  of  (2) ,  and  one  in  terms 

of  —P™  (l)  .  Thus,  the  Morillo  model  is  similar 

to  the  first  version  of  the  He  model  but  with 
both  m+n  odd  and  m+n  even. 

If  Eq.  (30)  is  substituted  into  Eq.  (17), 
the  set  of  ordinary  differential  equations  for  the 
velocity  coefficients  in  terms  of  the  pressure 
coefficients  for  skewed  flow  becomes, 

(35) 

In  axial  flow,  this  equation  becomes, 

[^]{»r}+[D']{«,"}=[D']{c}  <36> 

Equation  (35),  if  truncated  to  include  only  m+n 
odd  terms,  is  identical  to  the  form  of  the  He 
model  in  which  velocities  are  expanded  as 
p;(r).  Therefore,  the  He  model  is  implicitly 

included  in  this  new  model.  Furthermore,  since 
the  Pitt  model  is  a  special  case  of  the  He  model, 
the  original  Pitt-Peters  is  also  part  of  Eq.  (35). 

RESULTS 

In  this  section,  we  look  at  some  of  the 
results  from  the  newest  inflow  methodology  as 
compared  with  the  older  approach.  We  will 
consider  results  both  in  the  frequency  domain 
and  in  the  time  domain.  In  the  frequency 
domain  there  are  closed-form  numerical 
expressions  based  on  the  convolution  integral 
that  can  be  used  to  compare  the  solutions  with  a 
known  result.  There  are  also  closed-form 
solutions  in  the  time  domain  to  a  step  input. 
With  these,  we  can  determine  the  accuracy  and 
convergence  properties  of  the  methods. 

Figure  3  shows  results  in  the  plane  of  the 
rotor  disk,  through  the  y=0  center-line  both  on 
the  disk  (-l<x  <+l)  and  up  to  one  radial  length 
off  the  disk  (-2<x<+2).  The  normal  component 
of  flow  is  plotted  for  an  elliptical  pressure 
distribution  at  zero  frequency.  The  wake  skew 
angle  is  zero,  corresponding  to  purely  axial 
flow  as  in  hover  or  climb.  The  “exact” 


convolution  result  is  compared  with  both  the 
newer  Galerkin  approach  and  the  older  He 
model.  For  this  case,  all  three  results  are 
identical.  Only  the  real  portion  is  present 
because  of  the  zero  frequency.  Figure  4  shows 
the  same  results  but  at  a  reduced  frequency  of 
2.3.  Here,  the  new  approach  and  the 
convolution  have  identical  real  and  imaginary 
parts,  but  the  older  He  model  shows  some  error 
(of  the  order  of  10%).  This  is  the  worst  case  for 
discrepancy  between  the  old  and  new  models. 
The  difference  is  due  to  the  coupling  with  off- 
disk  inflow  modes  that  are  neglected  in  the  He 
approach. 

Figure  5  shows  the  same  result,  but  for  a 
cyclic  pressure  input  and  at  a  wake  sew  angle 
of  45°  and  zero  frequency.  Once  again,  the 
convolution  approach,  the  He  model,  and  the 
new  model  are  identical.  Figure  6  takes  the 
same  case  but  adds  a  reduced  frequency  of  4.0. 
Some  error  begins  to  show  up  in  the  new  model 
downstream  since  convergence  is  slower  for 
that  case,  and  only  ten  shape  functions  are 
included  in  these  results.  The  He  method  gives 
results  only  on  the  disk  (not  off);  but  it,  too 
starts  to  deteriorate  in  accuracy  toward  the  rear 
of  the  disk. 

Figure  7  gives  a  more  detailed 
examination  of  convergence  with  number  of 
even  terms  added  to  the  old  He  model.  The 
graph  gives  an  error  norm  (based  on  the  square 
of  velocity  errors  either  on  the  rotor  disk, 
-I<x<+1,  or  off-disk,  -2<x<+2)  for  the  normal 
component  of  velocity  in  axial  flow.  Ten  terms 
are  included  in  the  old  He  shape  functions 
(n+m  odd),  and  the  number  of  new  terms  (n+m 
even)  is  varied  from  zero  to  twelve.  (Note, 
since  terms  come  in  alternating  values  of  odd  or 
even  subscripts,  there  are  half  as  many  odd  or 
even  terms  as  the  value  of  the  highest 
subscript.)  One  can  see  how  the  error  quickly 
converges  to  zero  both  on  disk  and  off  disk. 
However,  if  too  many  even  terms  or  added, 
numerical  ill-conditioning  can  cause  the  error  to 
climb  back  up.  Therefore,  an  optimum 
approach  is  to  take  only  3/4  as  many  even  terms 
as  there  are  odd  terms. 

Figure  8  demonstrates  convergence  for 
flow  off  the  disk.  The  skew  angle  is  a 
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relatively  steep  75°.  Flow  is  plotted  one  rotor 
radius  above  the  disk  along  the  disk  centerline 
both  on  and  off  of  the  disk.  The  normal 
component  of  velocity  is  computed  for  10  odd 
and  10  even  terms  and  also  for  the  optimized  10 
odd  and  only  7  even  terms.  One  can  see  that 
there  is  improved  convergence  when  fewer 
m+n  even  terms  are  included  than  there  are 
m+n  odd  terms.  The  reason  for  this  is  that, 
while  the  m+n  odd  terms  (the  old  He  theory) 
are  very-well  conditioned,  the  m+n  even  terms 
are  poorly  conditioned.  Thus,  round-off  errors 
build  more  rapidly  in  these  new  terms. 
Fortunately,  the  on-disk  terms  give  a  good 
on-disk  result,  so  that  only  a  few  off-disk  are 
needed.  Then  flow  both  on  the  disk  and  off  the 
disk  is  satisfactory. 

Similar  results  are  obtained  for  the  other 
two  velocity  components  (azimuthal  and 
radial),  as  shown  in  Fig.  9,  which  gives  all  three 
components  of  flow  for  a  cyclic  input  at  zero 
degrees  skew  angle.  Peters-He  results  are  only 
shown  for  the  z  component  since  that  is  the 
only  component  for  which  the  model  can 
produce  results. 

In  addition  to  the  above  frequency 
domain  results,  correlations  have  also  been 
done  in  the  time  domain.  Figures  10  and  11 
show  results  in  the  time  domain  for  a  step  input 
in  pressure.  The  normal  component  of  flow  is 
plotted  for  the  case  of  axial  flow  and  for  a 
reduced  time  up  to  t  =7.  The  exact  solution  is 
compared  to  the  finite-state  results  with  10  even 
and  10  odd  terms.  Flow  is  plotted  at  the  disk 
center,  one-half  radius  out  from  the  center  and 
one  radius  out  from  the  center.  Virtually  exact 
agreement  is  obtained.  Figure  10  compares  this 
result  with  that  of  the  old  Peters-He  model  at 
the  rotor  center  and  at  the  rotor  edge.  The 
Peters-He  model  does  very  well  for  this  case. 

Non-Zero  Mass  Flux 

The  next  step  is  to  try  to  include  terms  with 
a  non-zero  net  mass  flux  along  any  blade,  d>” . 

This  is  based  on  the  fact  that  these  terms  yield 
infinite  kinetic  energy  in  the  flow  field.  Without 
losing  generality,  assume  a  pressure  distribution 


P  =  .  The  vertical  component  of  induced 

velocity  yields 


vt=4>S=P0»&°(‘'7) 


(37) 


and  therefore 


V,  = 


n 


on-disc 

-1  off-disc  <‘38') 

UJ 


For  off-disc  area,  if  Taylor  series  is  used  for 
large  r] 


then  the  kinetic  energy  per  unit  mass  of  the 
induced  flow  crossing  the  off-disc  rotor  plane 
(v  =  0,  a<jj<°° ) is  that 


Equation  (40)  fully  expresses  the  major 
discrepancy  of  including  non-zero  net  mass  flux 
terms  in  the  model:  Involvement  of  net  mass 
flux  terms,  which  happens  in  tremendous 
practical  cases,  will  theoretically  introduce 
infinite  energy  in  the  flow  field.  However,  we 
know  that  mass  source  rotors  exist  in  reality, 
and  that  they  do  not  have  infinite  kinetic 
energy.  This  is  due  to  two  facts:  the  media 
field  of  the  theory  does  not  dissipate  energy  at 
all,  and  the  whole  system,  including  the  media 
field  to  a  distance  of  infinity,  is  in  steady  state. 
In  reality,  first  of  all,  damping  of  air,  even 
though  very  small  under  the  condition  of 
helicopter  operating,  will  actually  dissipate  the 
energy  and  the  induced  velocity  decreases 
much  faster  away  from  the  rotor  than  those 
shown  in  Eq.  (39).  And  secondly,  the  operating 
period  of  helicopter  is  not  infinite.  Furthermore, 
the  study  of  induced  flow  is  not  the  final  goal  in 
the  study  of  helicopter’s  performance  -  the 
ultimate  goal  is  to  see  how  it  affects  the 
behavior  of  helicopter  blades,  fuselage  or 
personnel  on  the  ground  if  it  is  close  to  the 
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ground.  Therefore,  the  most  important  concern 
is  the  area  on  or  close  to  the  actuator  disc. 
Based  on  this  aspect,  if  the  results  have  good 
correlation  at  on-disc,  or  close-to-disc,  area,  the 
involvement  of  net  mass  flux  could  be 
tolerated. 

On  the  other  hand,  the  resultant  infinities  in 
[mc]  and  the  influence  coefficient  matrix  [£c] 

have  no  effect  on  the  damping  matrix  [dc]  , 

which  is  the  coefficient  matrix  of  the  pressure 
coefficient  vector  on  the  right  hand  side. 
Therefore,  it  might  be  possible  to  include  m-n 
in  r„m  but  not  in  the  velocity  potentials. 

To  include  net  mass  flux  terms,  i.e.,  P”  (v) 
and  Q"  (irj) ,  in  the  pressure  distribution,  the 
damping  coefficient  matrix  [dc]  on  the  right 

hand  side  will  be  required  to  have  extra 
columns.  Based  on  the  goal  that  the  new  model 
should  reduce  to  Eq.  (17)  if  net  mass  flux  terms 
are  not  considered,  these  new  entrees  are 
desired  to  follow  the  same  formulation  of  [dc]  . 

Thus,  we  use  the  same  formulas  as  in  the 
Appendix  but  with  m=n  allowed  on  the  right- 
hand  side. 

Considering  that  the  matrices  on  the  left 
hand  side  have  to  be  square,  the  extended 
momentum  equation  is  expressed  as 

[m'  ]{<C }+  [d‘  ][£'  ]"'  [m‘  ]{< } = [B‘  ]{C }  (41) 

where  [dc]  is  the  extended  damping  matrix. 

Results  with  Mass  Flux 

Figures  12-16  illustrate  the  on-disc 
optimized  frequency  responses  of  pressure 
distribution  P  =  with  various  system 
configurations.  In  all  plots,  circles  are  results 
from  the  convolution  integral;  triangles  are 
results  from  closed-form  solutions,  which  are 
available  in  limited  cases;  and  dots  represent 
results  of  the  proposed  state-space  model  by 
Galerkin  approach. 

As  discussed,  the  velocity  expansion 
excludes  any  terms  with  m  =  n ,  in  which 
m  —  n  —  0  is  an  extreme  case.  Notice  if  the 


pressure  distribution  is  P-  — O®  *  From  Eq.  (16), 
the  closed-form  solution  of  the  z-component  of 
on-disc  velocity  perturbation  is  uniform. 
However,  there  is  no  function  in  the  velocity 
expansion  that  has  property  of  uniformity  on- 
disc.  In  a  general  sense,  it  requires  a  large 
number  of  terms  in  the  velocity  expansion  to 
yield  a  good  approximation.  On  the  off-disc 
area,  the  Galerkin  approach  decays  into  the  far 
field  faster  than  the  exact  solution.  This  is 
expected  since  we  have  no  l/rj  terms  in  the 
velocity  expansions. 

Figure  13  is  the  frequency  response  for 
infinite  frequency.  Based  on  the  fact  that  the 
response  actually  becomes  zero  in  such  cases, 
ojvi ,  i  =  z,r,\jf ,  is  plotted  instead  of  v.  It  shows 
that  the  on-disc  and  far  off-disc  areas  have  good 
correlation  with  the  exact  solution,  yet 
relatively  large  oscillations  appear  at  the  edge 
of  actuator.  This  is  a  consequence  of  using 
continuous  expansion  functions  approximate  a 
continuous  velocity  expansion. 

Based  on  individual  error  analyses, 
optimized  results  of  frequency  response  under 
various  conditions  are  obtained  and  plotted  in 
Fig.  14. 

Figure  14  shows  response  of  a  higher 
frequency,  skewed-angle  flight.  Fig.  14  shows 
good  correlations  in  all  four  non-zero 
components.  This  is  a  direct  result  of  increase 
of  excitation  frequency. 

For  purposes  of  illustration  of  generality  of 
the  model,  Figs.  15  and  16  are  presented  to 
show  the  responses  of  pressure  distribution 
p  =  -<d|  with  zero  and  infinite  frequencies. 

From  all  the  plots,  it  can  be  seen  that  the 
radial  components  always  have  good 
correlation  with  the  exact  solution,  even  when 
the  axial  components  yields  large  error.  In  cases 
of  on-disc  optimization,  the  low  power  radial 
function  will  decrease  off-disc  errors  and 
therefore  yields  better  correlations. 

Thus,  the  theory  with  the  extended  [D] 
matrix  allows  for  non-zero  net  mass  flow. 

Effect  of  Wake  Curvature 

The  Pitt-Peters  and  Peters-He  models 
involve  the  wake  skew  angle  in  the  equations  of 
motion.  The  skew  angle  appears  in  the  form  of 
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X  =  tan(x/2)  where  %  is  the  skew  angle. 
Interestingly,  in  the  coupling  between  the  m-th 
and  r-th  harmonic,  X  appears  only  to  the 
powers  (m+r)  and  (m-r).  It  is  envisioned  that 
the  wake  curvature  angle  will  appear  in  a 
similar  way.  However,  since  the  wake 
curvature  is  small,  only  terms  to  the  first  power 
in  K  are  used.  These  have  been  determined  to 
involve  only  the  cross-coupling  between  the 
zeroeth  and  first  harmonics,  and  this  has  been 
done  for  the  He  model  but  not,  as  yet,  for  the 
newer  Morillo  model. 


Figures 


Fig.l  Coordinate  system. 


Fig.2  Volume  and  area  of  integration. 
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Fig.  3  Zero  frequency  response,  P  =  <P? ,  %=0°. 


Fig.  4  Frequency  response  (0=2.3,  P  =  <P° , 
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Fig.  5  Zero  frequency  response,  P  =  <p\,  X=45° . 
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x  :  Exact  Solution 

_ :  Galerkin  Approach  @  x=0 

- :  Galerkin  Approach  @  x=0.5 

Galerkin  Approach  @  x=l 


Fig.  10  Time  response,  P  =  &f ,  %=(f,  z=0,y=0. 


Fig.  11  Time  response,  P  =  ,  %=(f,  z=0,y=0. 
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Figure  12.  Frequency  response  of  pressure  distribution  P  =  -<I>°  with  co-  0.  All  other  components  of 
induced  velocities  are  zeroes.  Both  number  of  odd  and  even  terms  included  in  state-space  model 
are  11.  Highest  power  of  radial  polynomials  is  22.  Responses  are  evaluated  with  z  =  0“  (on  the 
actuator  plane),  j  =  0  (axial  inflow),  ^  =  0\180°  (along  fore-and-aft  axis).  Red  circles  present 
values  of  convolution  integral  results  at  locations;  blue  triangles  present  results  of  closed-form 
solution;  and  black  dots  present  results  obtained  by  Galerkin  approach. 


Figure  13.  Frequency  responses  of  pressure  distribution  P  =  with  0)  =  .  All  other  components  of 

induced  velocities  are  zeroes.  Number  of  odd  terms  included  in  the  state-space  model  is  10; 
Number  of  even  terms  is  3.  Highest  power  of  radial  polynomials  is  20.  Responses  are  evaluated 
with  z  =  0“  (on  the  actuator  plane),  z  =  0  (axial  inflow),  ^  =  0\180a  (along  fore-and-aft  axis). 
Blue  triangles  present  results  of  closed-form  solution;  and  black  dots  present  results  obtained  by 
Galerkin  approach. 
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Figure  14.  Frequency  responses  of  pressure  distribution  P  =  -<t>“  with  skew  angle  X  =  45°>  frequency 
co  =  4.  Number  of  odd  terms  is  12;  number  of  even  terms  is  4.  Highest  power  of  radial 
polynomials  is  Nr  =  24 .  Evaluation  is  performed  on  the  actuator  plane,  z  =  0' .  Plots  (a)  and  (b) 

are  real  and  imaginary  parts,  respectively,  of  vertical  component  of  induced  velocity.  Plots  (c) 
and  (d)  are  those  of  radial  component  of  induced  velocity.  Red  circles  present  values  of 
convolution  integral  results  at  locations;  and  black  dots  present  results  obtained  by  Galerkin 
approach. 
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Figure  15.  Frequency  response  of  pressure  distribution  P  =  -<!>!  with  co  =  0.  All  other  components  of 
induced  velocities  are  zeroes.  Numbers  of  odd  and  even  terms  included  in  state-space  model  are 
11  and  10,  respectively.  Highest  power  of  radial  polynomial  is  21.  Responses  are  evaluated  with 
z  =  0~  (on  the  actuator  plane),  (axial  inflow),  y/  =  0’ ,1 80s  (along  fore-and-aft  axis).  Red 

circles  present  values  of  convolution  integral  results  at  locations;  blue  triangles  present  results  of 
closed-form  solution;  and  black  dots  present  results  obtained  by  Galerkin  approach. 


Figure  16.  Frequency  response  of  pressure  distribution  P  =  -<t|  with  <y  =  «>.  All  other  components  of 
induced  velocities  are  zeroes.  Numbers  of  odd  and  even  terms  included  in  state-space  model  are 
12  and  4,  respectively.  Highest  power  of  radial  polynomial  is  23.  Responses  are  evaluated  with 
z  =  0~  (on  the  actuator  plane),  x~ 0  (axial  inflow),  (7  =  0°,  180"  (along  fore-and-aft  axis).  Red 
circles  present  values  of  convolution  integral  results  at  locations;  blue  triangles  present  results  of 
closed-form  solution;  and  black  dots  present  results  obtained  by  Galerkin  approach. 
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Appendix  I 

FORMS  OF  EXACT  SOLUTIONS 


9  v  9  v  _  9  P 

dt  9 z  9 z 

p(z,t)  =  f(t)p(z) 
p(z,t)  =  0  t<  0 


v(z»0)  =  v0(z) 


t>0 


Homogenous:  v(x,t)  =  v0(z-r) 

Step  Input: 

Outside  wake :  vow(z,t)=  P{z-t)~ p{z) 
Inside  wake:  vm  (z,  t ) = 2 P0  +  P(z-t)~  P(z) 


p(o_) 

p(P  )lK4KE 

z 

_ 

,  l 

General  Solutions :  p(o+  ) s  pn 

v„  =  v„  +  2/0P,  +  2P„  £7  f(r)dr  o  <z<t 


Add  Closed-Form  Integrals: 

vow  =f0P{z-t)-f{t)p(z)+  /(r)P(z-r  +  r)Jr 

V/w  =vovl, +2P0/(r-z)  0  <z<t 
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Vow  =-£*  f{t)P{z-t  +  t)dT 

v,w  =vow+2p(o+)f{t-z)  0  <z<t 

V/w=-=f/W^fc-^  +  ^+  Ho+)-P(oJ/(t-z) 

Note :  IW I  I  {'  /(r)  Pz  (z  ■ -  r  +  r]dr  =  f(t-z)  [p(o+  )  -  p( 0" )]/(*-*)+ =£  /  Wp<  (* " *  +  ^ 
Note :  =f  =  lim  [  +  f 

Jo  s — >0  Jo  J(-z+f 


Change  to  <f  =  z-t+r 


vo»  =-£,/(#-z+,)p«(#)‘if 

v„  =  -£  /(f-z  +  r)P,(#>if +/M) M«*)] 

— £  /(f-z  +  »)f>!(fV#  +  /((-z)  Ko-)+p(o*)] 
l  z_'  7  £  7 

- — y -  - Y 

continuous  added  mass  flow 


Outside  wake 

v  =  -  £  fin  ~  £  +  t)S/P{t])dT] 
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Check 

fXn-i+ifirWt) 

dt 

||  =  +  f  /'(>?-#+ t^P(>>)dri-f(i^P^) 

og  J-°° 

^  =  -/(<)Vf>(  f)  =  -V[/>(#,<)] 

df  d£ 

Check 

v(0) = -  £  f{r}~^P{v)dri 
T]<^=>f{neg)  =  0 

Inside  wake:  To  be  derived. 

V  =  ~  £  f{n  ~  £ + t)P{ri)  drj 

Appendix  II 

MOMENTUM  THEORY 


P0-Pa  =  \pvu{vu  +2vo) 


,v„  +  v0 
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Downstream:  P0  +AP+-j-pv%  =PD  +-^/o(v0  +vD)2 
_  Z _  £ _ 

P0-PD+AP  =  ^ pvD{vD  +  2v0) 


At  Rotor:  PDR  =  Pm  +  AP,  vUR=vDR=vR 


At  infinity: 


AP  =  2pvR(vR  +v0) 


At  Rotor: 

Pq-Pur  =|pvR(vy?+v0) 

P0  ~  Pdr  +  2 pvR  (vR  +  v0 )  =  (v 

-Po+Pdr  =|p^(2v0+3vr) 

Pdr-Pq  3vR  +  2vo 

ft- ft/fi  VR  2v0 

Hover: 

r D  rVR 

Low-Lift  Climb  v0  »  v 

Upstream: 

Pq-Pu  =Pvuv o 

Downstream: 

P0-PD+bP  =  pvDv0 

At  infinitv: 

AP  =  2pvRv0 

At  Rotor:  =  l 

P0  Pur 

Interactions 

The  Principal  Investigator  made  six  trips  to  Boeing  Helicopter  Company  in  Mesa, 
Arizona,  to  discuss  the  work  that  is  being  done  and  how  it  might  be  integrated  into  their 
simulation  codes.  The  presentations  were  also  made  to  the  Boeing  Foundation  in  St. 
Louis.  A  trip  was  also  made  to  Florida  Atlantic  University  to  discuss  the  inflow  model 
with  Professor  Gopal  Gaonkar.  Three  visits  were  made  to  AFTD  at  Ames  Research 
Center. 
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